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ABSTRACT 

We report on observations of recurrent jets by instruments onboard the Interface Region Imaging 
Spectrograph (IRIS), Solar Dynamics Observatory (SDO) and Hinode spacecrafts. Over a 4-hour 
period on July 21st 2013, recurrent coronal jets were observed to emanate from NOAA Active Region 
11793. EUV spectra probing plasma at transition region temperatures show evidence of oppositely 
directed flows with components reaching Doppler velocities of ±100 km s“^. Raster Doppler maps 
using a Si IV transition region line show all four jets to have helical motion of the same sense. 
Simultaneous observations of the region by SDO and Hinode show that the jets emanate from a source 
region comprising a pore embedded in the interior of a supergranule. The parasitic pore has opposite 
polarity flux compared to the surrounding network field. This leads to a spine-fan magnetic topology 
in the coronal field that is amenable to jet formation. Time-dependent data-driven simulations are 
used to investigate the underlying drivers for the jets. These numerical experiments show that the 
emergence of current-carrying magnetic field in the vicinity of the pore supplies the magnetic twist 
needed for recurrent helical jet formation. 

Subject headings: Sun: photosphere - Sun: chromosphere - Sun: transition region - Sun: corona - 
Sun: atmosphere - magnetic fields 


1. INTRODUCTION 


Since the disc overy of coronal jets in X-ray (e.g . Shi- 
bata et aL]|1992 ) and EUV (e.g. Chae et ar]|1999 ) imag- 


ing observations, there has been a growing body of ob 
servational and theoretical work investigating the phys 
ical mechanisms behind this phenomena. While the de 
tailed physical processes responsible for the acceleration 
of jet material depends on the local conditions (e.g. see 
Takasao et al. 2013), there is overwhelming evidence that 
magnetic reconnection is key for the impulsive energy re¬ 
lease associated with jets. As for the driving mechanism 
that allows for energy build-up, it has been reported that 
many jets are associated with emerging flux and/or flux 
cancellation eve nts in the photosphere. In the case of re- 
current jets (e.g.|Ghae et al.|1999|| Ghifor et al.|2QQ8l|Guo 

__ __ ' Chandra et aL||2015[ ) 

emanating from the same source region on the bun, an 
additional question is how the underlying driver leads to 
magnetic configurations that repeatedly erupt to produce 
jets of a homologous nature. 

We address the question posed above by performing a 
study of recurrent jets observed by multiple spaceborne 
observatories and by using data-driven simulations. The 
rest of the article is structured as follows. Section [2] 
presents observations of the jets in the transition region 
and corona by the Atmospheric Imaging Assembly (AIA, 
Lemen et al.|2Q12 Boerner et al.|[2Q12 ) onboard the So- 


cheung@lmsal.com 


lar Dynamics Observatory (SDO,|Pesnell et al.|2012) and 
by the Interface Re gion Imagmg Spectrograph (IRlS, |De 


Pontieu et al.|2Qi4 ). Sectionj^ presents photospheric ^ 


servations in and around the sourc e region of the jets by 
the Solar Optical Telescope (S OT, [Tsuneta et ^ j2QQ8 ) 
onboard t he H inode spacecraft ( [Kosugi et aL||2UU7[ r 
Section [4T] describes the evolution of the photospheric 
field in the region of interest as revealed by vector 
magnetogra ms from the Helioseismic and Magnetic Im¬ 
ager (HMI, jScherrer et al.||2012| |Schou et al. |2012|) on¬ 
board SDO. Sect ion |4. 2 1 presents results from simulations 
of coronal and chromospheric field evolution driven by 
HMI vector magnetograms. The physical implications of 
this study are discussed in section 


2. OBSERVATIONAL EVIDENCE FOR HOMOLOGOUS, 
HELICAL JETS 

The recurrent jets and the source region of these jets 
were simultaneously observed by IRIS, SDO and Hin¬ 
ode. Instruments onboard these satellites provide com¬ 
plementary coverage in wavelength, temperature and the 
spatiotemporal domains. Together they present an inte¬ 
grated picture of the magnetic and atmospheric environ¬ 
ment responsible for driving and initiating the set of re¬ 
current jets. The following sections discuss observations 
of the region of interest as seen by the various instru¬ 
ments. 


2.1. Photospheric and coronal observations by SDO 









































2 


Continuous full disk observations from SDO provide 
context about the environment in which the recurrent 
jets are generated. Fig. shows the large-scale structure 
of active region (AR) 11793 in the time period during 
which IRIS observed the recurrent jets. During this time 
range the AR is roughly centered at a Stonyhurst longi¬ 
tude and latitude of Wll deg and N23 deg, respectively. 

We refer to each of the four observed jets as Jl, J2, J3 
and J4 in chronological order. The four jets are shown 
in separate panels of Fig. Each panel consists of 
a HMI line-of-sight magnetogram (from the hmi.M_45s 
data series) overlaid with an EUV image from the AIA 
94 A channel. The aia_prep.pro routine in SolarSoft 
was used to align the full disk images from the two in¬ 
struments and to remap them to a common plate-scale 
of 0.6 arcsec per pixel. Inspection of the overlaid im¬ 
ages (and accompanying animation, available online) re¬ 
veal that the jets emanate from a strong plage region in 
the northeastern edge of the leading (negative) polarity 
patch of the AR. 

Though the jets are not identical, they possess strik¬ 
ingly similar features. Eirst of all, the ejecta in the jets 
are channeled into closed loops connecting the leading 
and trailing polarities of the AR. This is unlike some 
jets found in coronal hole regions, in which jet material 


is channeled into open magnetic field lines (e.g. Cirtain 


et al.||2QQ71 ISavcheva et al.||2QQ7[ |Moreno-Insertis et al. 
2QQ8[|PatsQiirakos et al.j 2008). Secondly, each of the jet 
events can be considered to possess a two-part structure, 
namely the inclined jet itself accompanied by the bright¬ 
ening of a compact c losed loop (or multiple closed loops) 


adjac ent to the jet (Shibata et al. 1994 Shimojo et al 


1996). It is worth noting that the jets can also be identi- 


hed m the other EUV channels. This likely implies that 
they have mult Thermal structure. 

In all four cases, the footpoints of the bright compact 
loops closest to the jets are found in the negative polarity 
network whereas the conjugate footpoints are located at 
a parasitic polarity patch within an adjacent supergran¬ 
ule. We will discuss the magnetic configuration of the 
environment around the jets in sections and 


2.2. Transition region observations by IRIS 

Between 11:34 and 16:34 UTC on 2013-07-21, IRIS ran 
a medium coarse 20-step raster observation program with 
150 repeats on the region of interest. The field-of-view 
of slit-jaw images (SJl) is 60" x 60" and the approximate 
pointing is indicated by the yellow box in the top right 
panel of Eig.[T] The raster step size is 2" so each spectral 
raster spans a EOV of 38" x 60". Each repeat of a spectral 
raster is accompanied by five slit-jaw images in each of 
the C II 1330, Si iv 1400 and Mg II k 2796 channels 
and one continuum image in 2832. Level 2 data was 
used for all of the following analysis. Slit-jaw images in 
the level 2 data product are dark-subtracted, flatfielded 
and geometrically corrected so that images from different 
channels are on a common grid. The same corrections 
are applied to NUV and EUV spectra. Eurthermore, the 
spectra are stacked as 3D raster cubes for convenient 
analysis. 

Eigure shows a Si IV 1400 slit-jaw image of Jl. 
The spectrograph samples along the vertical dark slit at 
X = 154 arcsec. EUV spectral line profiles are plotted for 


three different positions indicated by the blue, green and 
red cursors straddling the jet in the slit-jaw image. We 
first inspect the profiles of the Si IV lines at 1393.8 and 
1402.8 A. Both are transition regio n lines that form a t 
\ogT/K 4.9 (from CHIANTI 7.0, Landi et al. 2012). 
Eor this reason, the shapes of the line profiles are very 
similar. By visual inspection, it can be discerned that 
the Si IV line profiles at the three positions have very 
different bulk Doppler shifts. While the centroids of the 
blue and green profiles have blueshifts of ~ 80 and 30 
km s“^, respectively, the centroid of the red profile is 
clearly redshifted. The spectral profiles at all three loca¬ 
tions have wide wings with contributions beyond ±100 
km s“^ relative to the centroid positions. The spectral 
readout window of Si IV 1403 spectra includes the O IV 
1401 line, which forms at logT ^ 5.2. The signal from 
this line appears in the plotted Si IV 1403 spectra at 
Doppler shifts beyond —300 km s“^. 

The bottom panel of the right column in Eig[^ shows 
spectra for the C ll lines, which ar e expected to form at 
\ogT/K ~ 4.4 (from CHIANTI 7.0, |Laridi et al.|2012D . In 
the plots the rest wavelength of the red line (1335.71 A) 
was used to calculate the effective Doppler shift. The 
shapes of the C ll lines are much more complex than 
those of the Si IV lines. The profiles at the green and red 
cursor positions are bimodal. However the green profile 
has a higher amplitude peak on the blue side of the line 
while the red line has a higher amplitude peak on the 
red side. Even more complex is the blue profile, which 
has a trimodal shape. One may initially be tempted to 
interpret the bimodal profiles in terms of oppositely di¬ 
rected Doppler flows. However, we caution that the C ll 
lines are usually optically thick, so the dip in bimodal (or 
trimodal) line profiles likely results from opacity effects. 


nmodaij line pronies likely 
Fig-i is similar to Eig. 


but for J3. There are many 
qualitative similarities between the profiles in these jets. 
F^or instance the red and blue colored profiles here show 
bulk red- and blueshifts, respectively. In the case of J3, 
the lines are even broader (with contributions up to ±200 
km s~^ relative to the centroid position) and the pair of 
C II lines are now completely blended. The green profiles 
for the Si IV lines in Eig. show a bimodal structure. 
This is probably differe nt than the two-co mponent spec¬ 
tral profiles studied by Tian et al. (2012) in their study 
of EUV jets. In their case, they performed a double 
Gaussian fit to spectral profiles taken by the EUV Imag¬ 
ing Spectrograph (EIS) instrument onboard Hinode and 
found that one component corresponds to steady back¬ 
ground emission, while a blue-shifted second component 
is attributed to outflows from the jet. In our case (green 
Si IV 1394 profile in Eig. |^, one component has a bulk 
redshift while the other component has a bulk blueshift. 
Inspection of corresponding profiles at positions north 
and south of this location (i.e., as indicated by the blue 
and red cursors on the slit-jaw image) shows the pro¬ 
file to the north has a blueshifted component, while the 
profile to the south has a redshifted component. Since 
the profile sampled at the position of the green cursor 
is mid-way between the two, it is not surprising that 
the green spectral profile has both red- and blueshifted 
components. In contrast, the EUV line profiles in the 
post-impulsive phase (four minutes later, see Eig. are 
dimmer by almost two orders of magnitude, are much 
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Figure 1. HMI line-of-sight magnetograms (greyscale) overlaid with AIA 94 A channel images of the four recurrent jets. The ejecta from 
the jets are channeled into closed loops connected the leading and trailing polarities of AR 11793. The yellow box in the top right panel 
shows the field of view of IRIS slit-jaw images. 



Figure 2. Left: Si IV 1400 slit-jaw image of the first recurrent jet. Right: Spectral line profiles sampled at the blue, green and red cursor 
positions (all lying on the vertical slit) for Si IV 1403, Si IV 1394 and the pair of C ll lines at 1334 and 1336 (the rest wavelength of the 
latter is used as a reference). The spectra are plotted as functions of Doppler shift from their respective rest wavelengths. 
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Figure 3. Same as Fig. but for the third jet. The FUV profiles in this case are even broader, with lines having contributions from up 
to ±200 km s“^ Doppler ^ifts relative to the centroids. 
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Figure 4. Same as Fig. but 4 min later. The jet has already subsided and the FUV profiles are also much narrower than during the 
impulsive phase of the jet. 
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narrower and have zero mean Doppler shift. 

The FUV line profiles for different locations spanning 
the widths of J1 and J3 both indicate a bulk blueshift 
near the northern edge of the jet and a bulk redshift 
near its southern edge. Using raster scans, we investigate 
whether this is a pattern that pertains to all four jets. As 
per the above discussion, the Si IV 1394 line is simplest 
to analyze. Using this line, we calculate the zeroth and 
first moments 


Uine — dl 
('^line) — -^line 


( 1 ) 
( 2 ) 

where / is the spectrograph intensity (in DN sec“^ 
pixel” and v\os is the line-of-sight Doppler veloc¬ 
ity. /line is simply the integrated intensity across the 
line. Over the field-of-view of the IRIS rasters, the 
ratio / 1394//1403 is generally close to th e theoretical 
value of 1.95 as predicted using CHIANTI (|Landi et al. 
'2012|), which suggests the Si IV lines are optically thin. 
Adopting this assumption, we interpret ('i^une) to be an 
intensity-weighted mean Doppler velocity. 

Figure shows raster maps of /une and (vune) for the 
Si IV 139Tline for all four jets. While the amplitude of 
('^line) is larger in some jets than others, all four jets tend 
to have blue- and redshifts at their upper (northern) and 
lower edges (southern), respectively. This type of spatial 
pattern in Doppl er maps have previo usly been d etected 
in individual j ets ^|Pike fc Mason|19^ and su rges (Curdt 


hile Hike V 


& Tian 2011|). 

from the C oronal Dia gnostic Spectrometer ( CDS, [Har¬ 


rison et ar]|1995 ) and Curdt & Tian ( 2Q11|) used data 


from the Solar U I traviolet IVIeasuremen ts of Emitted Ra¬ 
diation (SUMER, [Wilhelm et al.|1995|) instrument, both 
studies used the same O V transition region line, which 
forms at \ogT/K ^ 5.4. In both studies, the spatial pat¬ 
tern of Doppler shifts was taken as evidence for helical 
motion. We adopt the same interpretation and take our 
IRIS Doppler maps as evidence for helical motion in all 
four of the recurrent jets. In a local Cartesian reference 
frame where I points along the jet (increasing height), 
the spatial pattern of {vnne) corresponds to rotational 
motion with vorticity cj/ = (V x v)i < 0. This implies 
the kinetic helicity of plasma motion associated with the 
jet is uiuji < 0 (by definition, ui is positive since it is the 
component of plasma flow along the jet direction). 

(5 IV lines observed in IRIS FUV spectra allow us 
to perform density diagnostics on the jet material. We 
used two different line ratios to measure the densities. 
The first is the ratio of the O IV 1401.1 and 1404.8 lines. 
The 1404.8 line is blended with a S IV line, so we used 
the S IV 1406.02 line and the assumption of optically 
thin emission to extract the intensity of O IV 1404.8. 
The O IV 1404.8 line is not always present in the spectral 
readout window of IRIS. However, in a number of slit 
positions where there is sufficient blueshift, we find the 
ratio O IV I40I.I/I404.8 to be in the range 4.0 — 4.5. 
For a temperature range of logT/K = 4.5 — 5.5 (derived 
from the ratio of S IV 1404.8 to S IV 1406.02), this ratio 
gives densities ranging from logUe/ cm”^ = 10.8 — II.0. 
Similarly, a ratio computed for the 1399 and 1401 lines 
of O IV has values in the range 0.29 — 0.35, which yields 


electron densities of logUe/ cm”^ = 10.8 — II.2. 

3. PHYSICAL DRIVERS OF RECURRENT HELICAL JETS 

Observations by SDO/AIA and IRIS establish the case 
that the jets are homologous and helical in nature. The 
four jets are homologous in the sense that they share 
substantial similarity in their observed spatial structure 
and evolution. This begs the question of the underlying 
driving mechanisms that lead to the initiation of the jets. 

There is a large body of work using numerical MHD 
simulations to study how emergence of ma gnetic flux 


V Shibata [19951 

Miyagoshi & Yokoyama 20041 Archon- 

tis et al. 

20051 G 

alsgaard et al.| 20051 Isobe et al.| 

2007 

Moreno- 

nsertis et al.||2008[ |Nishizuka et al.||2008||J 

Hegg- 

land et a 

.. 120091 Archontis et al.|20I01|Archontis & Hood 


2013 


1 

i 


ang et al.|2QI4[ ). Most of the simulations focus on 


single jets following reconnection between the emerging 
magnetic system with the ambient field. Recent work 
started to investigate how multiple jets can be emitted 
from the same source r egion. From a 3D MHD simula¬ 


tion of flux emergence, Archontis et al. (2010) reported 
that a series of reconnection events between the emerg¬ 
ing fl ux system and ambient coronal f ield le d to recurrent 


jets. Moreno-Insertis & Galsgaard (2013) performed a 


similar numerical experiment and toimd a succession of 
eruptions, some of which have physical properties that 
resembled the ‘standard’ type of jet while others were 
miniature flux rope ejections that that may be as sociated 
with so-called blowout jets (Moore et al. 2010). While 


( I998[ ) used dataT the Moreno-Insertis fc Galsgaard ~i |20I3 ) paper mentions 


that the erupting flux ropes in the simulation seem to 
rotate as if they were converting twist into writhe, the 
possible helical motion of the jets themselves were not 
studied. 


Fang et al. (2014) modeled solar jets by performing 
3D MHD simulations of twisted flux tubes emerging into 
a coronal layer with an ambient inclined field. They 
included magnetic field-aligned thermal conduction in 
their model, which provides the dominant mechanism 
for energy loss by plasma that has been heated to tran¬ 
sition region and coronal temperatures. They reported 
the existence of columnar jets consisting of plasma at 
a broad range of temperatures (chromospheric to coro¬ 
nal). Due to acceleration by the Lorentz force acting on 
reconnected field lines, the jet columns exhibit spinning 
motion. Synthetic intensity-weighted Dopplergrams (us¬ 
ing coronal lines) of a simulated jet column from a side 
perspective gives net Doppler shifts of ± 20 km s”^ on 
opposite si des of the column ar jet axis. 

Recently, 


ulation of t' 


Lee et al. (2015) performed a 3D MHD sim- 
le emergence of a strongly twisted flux tube 


from the convection zone into the atmosphere. They im¬ 
parted a density deficit distribution to the tube so that 
two segments of the tube would emerge and interact with 
each other as well as with the ambient inclined coronal 
field. The simulation yielded four episodes of twisted 
flux rope ejections carrying signicant mass loads away 
from the emerging flux region. These flux tubes untwist 
as they are ejected, resulting in signatures of torsional 
Alfven wave propagation. 

Apart from magnetic flux emergence, there is another 
means by which photospheric magnetic evolution can 
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Figure 5. Total intensity and mean Doppler velocity of the four jets as computed from IRIS observations of the Si IV 1394 line. In all 
four jets, there is a tendency for northern edge to be blueshifted while the southern edge is redshifted. This spatial pattern suggests all 
four jets are helical with the same (negative) sign of kinetic helicity. 


lead to jets, especially those with rotational motion. One 
possibility of generating h o mologous h elical jets was in¬ 


vestigated by Pariat et aL| (2009 2010). In their numer¬ 
ical MHD studies, the autnors considered the scenario in 
which a circularly shaped patch of magnetic flux is em¬ 
bedded in an environment with predominantly vertical 
field of the opposite sign. For convenience let us call the 
circularly shaped patch of magnetic flux the ‘parasitic’ 
pole. The initial potential field configuration has a coro¬ 
nal null with a fan that forms a quasi-separatrix (QSL). 
The QSL separates two magnetic volumes, one consist¬ 
ing of the set of closed field lines connecting the parasitic 
pole with its surroundings, and an exterior volume con¬ 


sisting of purely open magnetic field lines. About the 
axis of symmetry of this parasitic pole, they applied ro¬ 
tational motion at a fraction of the local Alfven speed to 
mimic twisting of the field due to horizontal photospheric 
flows. The numerical simulations showed that, given suf¬ 
ficient cumulative rotation (greater than one turn), the 
magnetic system is driven to a state where reconnection 
of the magnetic field allows it to impulsively release the 
stored magnetic energy. A consequence of the relaxation 
is the formation of a helical jet. Subsequently, continued 
application of rotational driving at the bottom boundary 
led to the formation of similar helical jets. 

The driving mechanism considered by Pariat et al. 
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(2009, 2010) is in principle different to that of emerging 
flux in that only horizontal motions at the photosphere 


the vicinity of the parasitic pore. These bright grain pairs 
are separated by a dark lane with a length of 1 — 2 Mm. 


also drive systematic horizontal flows that shear already 

flux (e.g. Strous & Zwaan 

1999 Cheung et al. 2008 

emerged field (Manchester IV|200ll 

Magara & Longcope 

Guglielmino et aL|20I0 ) and their presence in the vicinity 

20031 IManchester IV et aL||2004[ |1V 

lagarall'JUOtil 

Cheung 

of the parasitic pore suggests flux emergence may play a 

et al.|2010t|Fang et al.|2010| 2012 Cheung & Iso 

be|20141 

role in driving the recurrent jets. 


to injection of magnetic energy and magnetic helicity for 
driving eruptive phenomena, including jets and CMEs. 
So the physical mechanism driving the recurrent jets re¬ 
ported here may not be exclusively due to one scenario 
or the other (i.e. flux emergence vs. shearing/rotational 
motion). 

To investigate the underlying physical driver(s) that 
cause the observed helical jets, we examine observations 
of the lower atmosphere by Hinode/SOT and SDO/HMI. 

3.1. Lower atmospheric observations by Hinode/SOT 

Figure M shows SOT observations of the jet emitting 
region ofAR 11793. The left panel shows the continuum 
intensity image from a Hinode/SP raster scan. The cen¬ 
tral panel shows the vector magnetogram from a Milne- 
Eddington S tokes 


et al. 


inversion (|Lites & Ichimoto|2013| 
right panel shows a broadband 


2013). The right pan^ 


Lites 


inter 

image from the Ca ll H channel during ocurrence of the 
third jet. The latter shows a set of closed loops and a set 
of inclined loops pointed toward the northeast direction. 
The morphology o f the loops is som ewhat reminiscent of 
the jet studied by Liu et al. (2009| see figure 1 of their 
article). In their paper, they studied SOT observations 
of a jet using Ca ll H images alone. From the morpholog¬ 
ical evolution, they concluded that the apparent motion 
of the jet material was consistent with a helical jet with 
untwisting magnetic field lines. In the case they studied, 
the jet was observed off the solar limb so the two-part 
structure comprised of the closed loops and the inclined 
open loops were not contaminated by background emis¬ 
sion. In our case, contribution to the SOT Ca ll H chan¬ 
nel by emission in the upper photosphere makes it harder 
to delineate the specific features. Still, the morphologi- 
cal similarities with the case studied by Liu et al. (2009) 
support their conclusion that their jet exhibited helical 
motion. 

The eastern (left) ends of the closed loops in the Ca ll 
H image appear to be anchored at a compact patch of 
positive polarity field located at (x^y) = (165", 287"). 
Inspection of the accompanying continuum image shows 
this positive patch to be a pore. This positive polarity 
pore is embedded in the interior of a supergranule. In 
contrast, the network field surrounding this supergran¬ 
ule is predominantly negative (the same polarity as the 
nearby leading spot). In this sense the pore is a parasitic 
pole. Inspection of the transverse field (i.e. the ampli¬ 
tude of the plane-of-sky component \Bt\) shows the pres¬ 
ence of strong horizontal field {Bf > 700 G) on the west¬ 
ern (right) side of the pore. Furthermore there is a mod¬ 
erate transverse field {\Bt\ ^ 300 G) pervading a large 
fraction of the sugergranular cell. The left-right asym¬ 
metry of the Bt distribution about the pore is suggestive 
of electric currents associated with a non-potenti al m ag- 
netic field configuration (more on this in section 4.1). 

Inspection of the temporal sequence of Ga ll H images 
indicates the presence of bright grain pairs that form in 


4. MAGNETIC FIELD EVOLUTION 

In this section we investigate how the magnetic field in 
the source region of the jets evolve, and how this leads to 
recurrent jet production. For this purpose we use HMI 
vector magnetograms to drive simplified numerical mod¬ 
els of chromospheric and coronal field evolution. 

4.1. HMI Vector Magnetograms 

The Hinode/SP vector magnetogram shown in Fig. 
suggests the presence of electric currents in the photo- 
spheric field around the parasitic pore. Due to Hinode 
telemetry limitations, only one SP map is available in 
the interval containing the four jets. So for the purpose 
of inspecting the evolution of the photospheric magnetic 
field, we used vector magnetograms from HMI instead. 

Each HMI vector magnetogram is produced by a 
Milne-Eddington inversion of IQUV Stokes maps tempo¬ 
rally in terpolated over an apo dization window spanning 
1350 s ( Hoeksema et al.|[2Q14 ). This is done to increase 
the signal-to-noise ratio of the Stokes parameters and to 
filter out p-mode oscillations. While HMI vector mag¬ 
netograms are not instantaneous representations of the 
photospheric field at any given time, the pixels in the 
same magnetogram are co-temporal. HMI provides vec¬ 
tor magnetograms of the full AR at a regular cadence of 
one frame per 12 min. We use vector magnet ograms from 
the HMI data series hmi. sharp_cea_720s (Sun 2013). 
This series provides vector magnetograms in Ar patches, 
such that the magnetic field vectors have been disam¬ 
biguated, transformed and remapped onto a cylindrical 
equal area grid. The magnetic vector is expressed as 
{Br, Be, B(f)), corresponding to the radial, longitudinal 
and latitudinal components, respectively. Since our main 
region of interest is relatively small {L ^ 30 Mm) , we 
ignore the effects of curvature and use the following map¬ 
ping to a local Gartesian coordinate system Bq Bx, 
B(f) By and B^ ^ B^. This enables us to compute the 
vertical current density 


dB,, dBx 


dx 


dy 


( 3 ) 


Inspection of maps of jz during the jet-emitting period 
reveals a persistent patch of positive jz near the para¬ 
sitic (positive polarity) pore. Fig. shows HMI vector 
magnetograms in the neighborhood of the pore at 09:00, 
10:46, and 12:34 UT. The left column shows the verti¬ 
cal component Bz. The positive polarity pore is roughly 
centered at {x,y) = (—2,0) Mm in all three snaphots. 
The horizontal magnetic field can be decomposed into the 
sum of a potential component (calculated from Bz) and 
a current-carrying component. These are respectively 
shown in panels in the middle and right columns. A per¬ 
sistent patch of current-carrying field can be found on the 
west side of the pore at {x,y) = (—1,0). This feature is 
also found in the middle panel of Fig. which shows 
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Figure 6. Hinode SOT observations of the jet emitting region of AR 11793. Left: Continuum image from a Hinode/SP raster scan. Center: 
Vector magnetogram of the same field of view. The line-of-sight component {Bi) of the magnetic field is denoted such that blue and red 
denote positive and negative polarities, respectively. Green color coding shows the strength of the transverse component {Bt). Overplotted 
lines show the orientation of the transverse field. Right: Broadband Filter (BFI) image in the Ca II H channel during the occurrence of 
the third jet. 

the Hinode/SP vector magnetogram of the same region. 

The Hinode/SP map provides a more accurate field mea¬ 
surement due to higher spectral coverage and resolution. 

The strength of the current-carrying horizontal field in 
this patch reaches 1 kG and exceeds the field strength 
of the current-free counterpart. This indicates that the 
magnetic field on the west side of the pore is strongly 
twisted. Inspection of the Bz time sequence from HMI 
shows that magnetic flux is emerging in this area. The 
orientation of the emerging flux is such that negative and 
positive polarity field are migrating in roughly the north 
and south directions, respectively. So both Hinode/SP 
and SDO/HMI vector magnetograms give evidence for 
the emergence of twisted field in the vicinity of the pore. 

In the following section, we present numerical simulations 
to investigate how this photospheric driving is related to 
the phenomena of recurrent helical jets. 


4.2. 


Numerical Experiments using a Data-Driven 
Magneto frictional Model 


We use a time-dependent magnetofrictional (Yan; 


et al.||198^ [Craig fc Sneyd p^86| 

2QQQp model to carry out data-d: 


van Ballegooijen et a. 

to carry out data-driven numerical experi- 

ments of coronal field evolution. Under this model, the 
fluid velocity v appearing in the induction equation is as¬ 
sumed to be proportional to the local Lorentz force j x B. 
This leads to an evolution of any arbitrary magnetic con¬ 
figuration to relax toward a force-free field. Magnetofric¬ 
tion has been used to model the formation and evolution 


of filaments (Mackay et al.|2000 

van Ballegooijen[[2004i 

Mackay & van Hallegooiien|2UU6l 

20091 [Yeates et al. 2007, 

2008 Yeates & Mackay|2009p , the coronal held above the 

quiet Sun magnetic carpet (jlVLeyer et al. 

2012|) and the 

evolution of ARs (Cheung DeRosa [20 J 

Cibb et al. 


Following Cheung & DeRosa ([2012^, we use a Carte¬ 
sian magnetofriction code that solves for the vector po¬ 
tential H, namely 




where H = VxH, ji=VxH and 

v = fjxB. (5) 

The magnetofrictional coefficient n is given by 

v = voB\l-e-^/^), ( 6 ) 

w here z^n = 10 s Mm~^ and L = 1.7 Mm. As described 
in [Cheung fc DeRosa (2012), the code uses a staggered 
grid (Yee mesh) such that A and j are defined at cell 
edges, B is defined on cell faces and v is defined at cell 
centers. The code has been updated to use a van Leer 
slope limiter (van Leer 1977) to interpolate v onto cell 
edges when computing the —v x B electric field. We 
find that this scheme provides better numerical stabil¬ 
ity while being less diffusive than explicitly imposing an 
anomalous resistivity. 

4.2.1. Initial and boundary conditions 

The observed recurrent jets occur in the neighborhood 
of ambient inclined field as part of a set of AR loops that 
connect the leading and following polarities of AR 11793. 
To capture the large-scale magnetic connectivity, we use 
a computational domain sufficiently large to encompass 
the entire AR. The domain spans 248 and 131 Mm in the 
X (longitudinal) and y (latitudinal) directions, respec¬ 
tively. The bottom boundary is located at 2 ; = 0 and the 
top boundary is located at 105 Mm. The horizontal and 
vertical grid spacings are 364 and 547 km, respectively. 
The initial condition is a potential field of the AR com¬ 
puted for the magnetogram at 2013-07-21T06:12 UT (5 
h 22 m before the start of the IRIS observation) and the 
simulations are e volved forward in time fro m that state. 

As discussed in [Cheung & DeRosa[ ([2012[), the bottom 
boundary condition tor the magnetofrictional model is 
given by the transverse components of the photospheric 
electric field, namely —dfAx = Ex and —dtAy = Ey. 
The retrieval of the horizontal electric field from vector 


(4) magnetograms is a difficult inverse problem (Fisher et al 
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Figure 7. SDO/HMI vector magnetograms of the parasitic pore and its surroundings. Left: Vertical component of B (greyscale saturated 
at ±200 Mx cm“^). Middle: Horizontal components of current-free (i.e. potential) part of B. Right: Horizontal components of the 
current-carrying part of B. A strong, persistent patch of current-carry field is found on the west side of the pore. Comparison with the Bz 
distribution shows this current-carrying patch is coincident with an emerging fiux region (just northwest of the parasitic pore). Contours 
for Bz = 500 and 1000 Mx cm“^ are shown on all panels to indicate the position of the pore. 


2Q1Q[ |2Q12|). As recently demonstrated by Kazachenko 


et al.| (2014), the incorporation of Doppler hows as a 


constraint leads to inversion results that accurately re¬ 
produce electric fields (and the associated Poynting flux) 
in an anelastic MHD simulation. The validation of this 
method for use with vector magnetograms at the reso¬ 
lution and sensitivity of HMI on photospheric magnetic 
structures is work in progress. 

Instead of attempting to perform a faithful retrieval 
of the photospheric electric field for this problem, we 
use a different method to compute electric fields with a 
given assumption. Given the sequence of input vertical 
magnetograms its relation to = {E^^ Ey) is given 
by the vertical component of the induction equation 


dB^(x,y) 

dt 


= • (V X Eh ). 


( 7 ) 


In order to solve for E^^ another relation must be speci¬ 


fied. Specifying the horizontal divergence of the electric 
field: 

D{x,y) = Vh • Eh- ( 8 ) 

fully constrains the problem. In our numerical experi¬ 
ments, we assume 


D{x,y) = j^{x,y)Uo, (9) 

where jz = {\/ x B)z is the vertical current density com¬ 
puted at the photosphere. This choice of the functional 
form for D is motivated by the following scenario. Con¬ 
sider an axisymmetric twisted magnetic flux tube that is 
invariant along its axis and let the tube axis be parallel 
to the vertical direction z. Let the tube rise vertically 
with V = Uqz. It can be shown that the divergence of 
the —V X B electric field corresponding to this motion 
is given by Eq. By adopting Eq. ([^, we have cho¬ 
sen to inject twist via the bodily emergence of twisted 
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field (c.f. Leka et al. 1996). For driving the numerical 
experiments described here, jz is computed from HMI 
vector magnetograms and Uq is a free parameter (with 
dimensions of velocity). Varying Uq changes the effec¬ 
tive injection speed of magnetic twist (as described by 
jz) into the computational domain. 

Another possible way to inject twist is by means 
of shearing and rotational motions in the photospheric 
plane. In MHD models of the emergence of twisted flux 


exerted by the Lorentz force ( 

Longcope & Welsch 

2000: 

Manchester IV et al. |2004l 

Magara 

|2006[ |F'an| 


Cheung et al. 20101 |Fang et al. 2012 

Cheung & Isobe| 


current-carrying field and rotational motions in the pho- 
tospher e ar e likely present in our region of interest (see 


section 3.1). For simplicity, however, we assume in our 
numerical setup that twist injection is due to the former Q 

4.2.2. Results of Numerical Experiments 

We carried out three numerical experiments, one each 
for Vo = 0, Vo = 1.1 and Vo = 2.2 km s“^. For the 
run with Vo = 0, Eh is decoupled from the photospheric 
distribution of jz- In this case there is no systematic 
current injection and we do not find any magnetic field 
evolution resembling the helical motion of the observed 
jets. For Vo = 1.1 km s“^, we find repeated episodes of 
twisting followed by untwisting in the model magnetic 
field that shares qualitative similarities with the recur¬ 
rent jets. Fig.j^shows snapshots of two such homologous 
episodes of evolution from the model. Row (a) shows Bz 
at the bottom boundary of the domain, which is con¬ 
strained to match the Bz observed by HMI. Row (b) 
shows synthetic chromospheric “magnetograms” (i.e. Bz 
sampled at 2 : = 4.5 Mm). We chose to sample the mag¬ 
netic field at a height that is above the nominal height 
range associated with ID mo dels of the solar chromo 
sphere (^ 1 to 2 Mm , e.g. see Withbroe & Noyes 1977 


Vernazza et al. 1981) since 3D radiative MHD simula¬ 


tions show that magnetic flux emerging into the atmo¬ 
sphere can lift the ch romosphere and transition re gion 
layers by a few Mm (Martmez-Sykora et al. 2008|). In 
both cases, the earliest chromospheric magnetogram (i.e. 
panels (bl) and (b4) ) shows a positive polarity feature 
pressed against a negative polarity feature in a ying-yang 
pattern. Inspection of the horizontal vectors and the dis¬ 
tribution of jz at that height shows that the field near 
the polarity inversion line is sheared and carries current. 
When we inspect subsequent vector magnetograms, we 
find the horizontal vectors have, on average, rotated in 
a clockwise direction. This is a result of the untwisting 
of the magnetic field by a clockwise rotational flow. The 
vertical component of vorticity for such a flow is nega¬ 
tive, consistent with the sign of vorticity of the observed 
recurrent jets. 


^ For the specific case of an axisymmetric flux tube with az¬ 
imuthal held Bgir) = qrBi(r), where q is the twist parameter 
and Bl (r) is the longitudinal component of the magnetic held, one 
can pick parameters for the two scenarios (i.e. bodily transport of 
twisted held and rotational motion) that result in identical bound¬ 
ary conditions. For such a tube rotating about its axis with an 
angular velocity cao, it can be shown that V • Eh = —q~^(jOojz, 
which is equivalent to Eq. ^ with Uq = —qujQ. 


To visualize the magnetic field lines, we calculate a 
scalar proxy emissivity e{x, y, z) using the following pro¬ 
cedure. For each field line traced from a position at the 
photospheric {z = 0) boundary, we compute the follow¬ 
ing field-line averaged quantity: 

(WD) = L~^ [ FLorentz ' vdl, (10) 

Jo 

where L is the length of a field line and Fporentz * v rep¬ 
resents the rate of work done by the Lorentz force. In 
a full MHD model, Fporentz * ^ appears as a source term 
in the kinetic energy equation and as a sink term (with 
negative sign) in the magnetic energy equation. In a 
magnetofrictional model, 

FLorentz ' V = [(V X F) X F]^ > 0. (11) 

47rz^ 

so that (WD) >0. If a field line crosses one of the side 
or top boundaries of the computational domain, we set 
(WD) = 0 so that the field line is not emissive. A mag¬ 
netic field line will cross a number of cell elements in the 
computational domain. For each of these cell elements, 
we increment the local value of the emissivity by 


de = G{WD)AxAy, 


( 12 ) 


where G is some arbitrary scale-factor (here we use 
G = 1). Row (c) shows integrals of the resulting proxy 
emissivity along vertical lines-of-sight. This method al¬ 
lows us to inspect a large number of fields lines (here we 
traced 4 field lines per pixel), each lit up according to 
the field-line averaged rate of magnetic energy loss due 
to work done by the Lorentz force. The images shown 
in row (c) are log-scaled so they do not depend on the 
specific value of the constant G. 

Fig.i shows a 3D visualization of the magnetic struc¬ 
ture from the same data-driven simulation {Uq = 1.1 
km s“^) at 11:02 UT. The topological structure in this 
case is similar to t hose in scenarios examined by [Pariat 
et al. (2009, 2010). In their numerical experiments, a 
parasitic polarity is embedded in an ambient fie ld that 
is either ve rtical (Pariat et al. 2009) or inclined (Pariat 
et al.| 2010). In both cases, the model coronal field has 
a null point. Associated with the null point is a spine 
field line connecting the parasitic polarity with the null 
and a set of magnetic field lines for ming a fan locus (i.e . 
a so-call fan-spine_topology, see e.g. Parnell et al.||1996). 
As shown in Fig. the same type of magnetic skeleton 
is found in the data-driven model. In this case, the am¬ 
bient field (with polarity opposite to that of the pore) 
is concentrated at the supergranular boundary, so the 
fan field lines emanating from the coronal null are con¬ 
nected to network flux. This type of topology is similar 
to the one inferred fo r th e (un)twisting jet event studied 
by Guo et al. (2013) and Schmieder et al. (2013). That 
event also occurred above a parasitic patch embedded in 
a supergranule at the edge of an AR. 

As demonstrated by the nu merical MHD experiments 
by Pariat et al. (2009 2010), persistent rotation of the 
parasitic polarity leads to a series of homologous helical 
jets. The simplistic MF model used in our data-driven 
simulations do not solve for the continuity, momentum 
nor energy equations so the model does not yield plasma 
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Figure 8. Two examples of untwisting magnetic field in the data-driven model with Uq = 1.1 km s ^. Columns (1) - (3) show a temporal 
sequence of one such episode while columns (4)-(6) show a homologous episode later in time. In this local coordinate system the parasitic 
pore is centered at the origin. Row (a) shows the photospheric vertical field Bz{z = 0) scaled between ±800 Mx cm“^. Row (b) shows 
synthetic chromospheric vector magnetograms at a constant height of 2 : = 4.5 Mm. Greyscale shows Bz scaled between ±100 Mx cm“^ 
and red arrows show the horizontal field at the same time. In the left and right sets of panels, the green arrows indicate the horizontal field 
at 10:48 and 12:45 UT, respectively (i.e. they show the horizontal field at the beginning of each set). The clockwise rotation of the arrows 
in both cases shows the magne tic fi eld evolved with a clockwise rotation. Row (c) shows a visualization of the field lines according to the 
proxy emissivity model of Eq. 


ejections. However one can still examine the magnetic 
evolution in the data-driven model and find similarities 
between magnetofrictional evolution and MHD evolu¬ 
tion. Fig. shows two snapshots from a 3D visual¬ 
ization of the magnetic evolution in magnetofrictional 
model with Uo = 1.1 km s“^. At 12:55 UT in the model, 
the pink field lines reveal a twisted flux rope connecting 
the parasitic pore with the northwestern edge of super¬ 
granule boundary. At 13:37 UT, the field has evolved so 
that the pink field lines (traced from the same positions) 
trace out a twisted bundle of inclined field. The sign 
of magnetic twist in this flux bundle is consistent with 
the sign of rotation revealed by IRIS Doppler mean shift 


maps shown in Fig, 
bundle results in a 


0 

i L( 


That is, the sense of twist in this 
(Orentz force that drives rotational 


motion with cj/ < 0 (see sec tion 


_ (s( _ 

In thejPariat et al^(f2009|) model, the parasitic polarity 
was embedded in an ambient field that is purely vertical. 
For this setup, they found that the critical number of 
windings needed to be injected in to the system to form a 
helical jet to be = 1.4 (see also Rachmeler et al.|2QlQ ). 
When the number of windings injected by surface rota- 
tion reached this value, the system underwent a helical 
kink instability, which broke the azimuthal symmetry of 
the system and generated a helical jet. The exact value 
of N depends on the geometry of the system. When the 
ambient open field is inclined, the azimuthal symmetry is 
no longer present and the critical threshold for injection 


2 . 2 ). 


of twist is lower, with a value of = 0.85 ±0.1 (jPariat 

eF^pMot - 


The average length of time required for twisting the 
field before jet-like reconfigurations occur in th e data- 
driven mod els is consistent with the findings of [Pariat 


et al. (2010). To estimate the windings injected per unit 
time m the data-driven simulations, consider an axisym- 
metric twisted flux tube with longitudinal and transverse 
components given by 


= (13) 

Bt{r) = ^Boe-^"/^\ (14) 


Here r is the radial distance from the axis of symmetry, 
A is the dimensionless twist parameter and R is the char¬ 
acteristic radius of the flux tube. This magnetic config¬ 
uration consists of helical field lines that form concentric 
flux surfaces about the tube axis. The field lines have 
magnetic pitch such that the number of turns about the 
axis (A) over an axial distance d is independent of radial 
distance r and is given by: 


N = d 


A 

2 ^' 


(15) 


Consider the scenario in which magnetic twist is injected 
into the corona by a vertically aligned twisted flux tube 
rising through the photospheric layer {z = 0) with speed 
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(a) Bz dit z = 0 Mm 




(a) Close-up view of magnetic flux rope at 12:55 UT 


(b) Bz dit z = A Mm 



(c) at 2 : = 8 Mm 



(b) Magnetic fleld at 12:59 UT 


(c) Magnetic fleld at 13:37 UT 



Figure 9. The magnetic conflguration in the neighborhood of the 
parasitic pore in the data-driven model (with Uq = 1.1 km s“^) at 
time 11:02 UT. Magnetic fleld lines traced from seed points in the 
vicinity of a coronal null point reveals a fan-spine topology, with 
the spine connecting the null point with the underlying parasitic 
polarity. 


Uq. Over a time period At, the number of turns injected 
is N = UoAtX{27TR^^ . By inspecting the HMI vector 
magnetograms (fig. 1^, we find that the horizontal com¬ 
ponents of the current-carrying field in the vicinity of 
the parasitic pore to be comparable to the vertical field 
strength inside the pore (approaching 1 kG). So for rough 
estimates we can take A ^ 1. Taking the size of the pore 
8iS R ^ 1500 km, we find that the time taken to reach 
the critical threshold of N = 0.85 to be At = 2.2Uq^ hr, 
where Uq is in units of km s“^. Within the 4-hour period 


Figure 10. Jet-like magnetic evolution in the data-driven model 
(with Uq = 1.1 km s“^). Pink magnetic fleld lines are traced from 
a stationary grid of points. At 12:55 UT, they show a twisted flux 
rope structure. At 13:37 UT, fleld lines traced from the same set 
of points reveal a set of inclined twisted fleld lines aligned with 
the background inclined fleld. Green fleld lines are traced from 
2 ; = 0 from the parasitic polarity. The semitransparent orange/red 
surfaces in panels (b) and (c) indicate regions of strong current 
density. 


in the simulations (10:00 UT and 14:00 UT), the number 
of jet-like reconfigurations in our data-driven simulations 
with Uo = 0, 1.1 and 2.2 km s~^ are 0, 2 and 4, respec¬ 
tively. So the amount of twisting required to drive jet-like 
episodes in the data-driven model is consistent with w hat 
is required in the MHD model of Pariat et al. (2010). 

There are limitations to the rnagnetolrictional rnod- 
els. Although they give us physical guidance about how 
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magnetic energy accumulates and how the magnetic field 
relaxes by unwinding, the amplitude of velocities in the 
models is dependent on the free magnetofrictional pa¬ 
rameter vq. For smaller values of the relaxation ve¬ 
locities have higher amplitude. The magnetofrictional 
models were not fine-tuned to reproduce the Doppler 
speeds in the observations (generally in our model the 
speeds are lower by a factor of a few). Furthermore, 

the relaxation velocity {v oc j x B) is always perpen¬ 
dicular to H, so the model cannot give predictions of 
the outflow speeds of jet material (from the IRIS slit- 
jaw image sequences, the outflow speeds of the jets reach 
beyond 100 km s“^). Another limitation is the lack of 
a treatment of thermodynamics quantities such as tem¬ 
perature and density. To overcome these limitations 
would require data-driven, fully-compressible MHD sim¬ 
ulations to be performed in future studies. Another 
valuable exercise would be to compare the amount of 
twist injection needed in MHD simul ations of flux emer- 


gence 

(those yielding helical jets, e.g. 

Archontis & Hood 

2013 

. 

IMoreno-Insertis & Galsgaard| 

2013J pang et al. 

2014 

iLee et al. 2015) with the resu. 

ts of Fariat et al. 

(|2U1U 



5. DISCUSSION 

Over a four-hour period on July 21st 2013, recurrent 
jets emanating from NOAA AR 11793 were simultane¬ 
ously observed by IRIS, SDO and Hinode. Doppler shift 
maps in the IRIS Si IV 1394 A transition region line 
shows all four jets exhibiting helical motion of the same 
sign. The IRIS Doppler shift maps share considerable 
resemb lance to synthetic Doppler maps in Fang et al. 
(2014), who carried out MHD simulations of jets result¬ 
ing from the interaction of a twisting flux tube emerging 
from the solar convection zone into a coronal layer with 
ambient inclined field. 

Photospheric vector magnetograms from Hinode/SOT 
and SDO/HMI show that the source region of the ho¬ 
mologous jets consists of predominantly negative polar¬ 
ity field concentrated at the boundary of a supergran¬ 
ule. Embedded inside the supergranule is a parasitic pore 
with positive magnetic flux. This type of photospheric 
flux distribution gives a coronal field with a spine-fan 
topology, wh ich is common in 3D MHD models of coro¬ 
nal jets (e.g. [Moreno-Insertis et al.||2QQ8 ). 

Photospheric vector magnetograms from Hinode/SP 
and SDO/HMI show a persistent current-carrying mag¬ 
netic configuration in the vicinity of the parasitic pore. 
Furthermore, the temporal sequence of vector magne¬ 
tograms from SDO/HMI shows evidence for the emer¬ 
gence o f ma gnetic flux in this current-carrying region (see 
section 4.1). 


To investigate the driving mechanism for the homolo¬ 
gous helical jets, we performed a number of data-driven 
numerical simulations. All of the numerical simulations 
are driven by a bottom boundary condition that matches 
the evolution of the photospheric-as observed by HMI. 
The temporal sequence of Bz indicates flux emergence is 
in progress during the time period when helical jets are 
observed. The occurence of helical jet-like evolution in 
some simulation runs and not others implies (for this 
particular case) that the increase in unsigned magnetic 
flux {\Bz\) associated with emerging flux is not a suf¬ 


ficient condition for helical jet formation. What vector 
magnetograms (Fig. reveal is that the emerging flux 
is current-carrying (i.e. has magnetic twist). In the nu¬ 
merical experiment for which the driving electric fields at 
the bottom boundary are completely decoupled from the 
photospheric jz distribution, we find no helical, jet-like 
reconfigurations in the magnetic field model. In cases 
where twisting is imposed (i.e. the driving electric fields 
are coupled to jz)^ the number of jet-like episodes within 
a 4-hour period increases linearly with the injection pa¬ 
rameter Uq (see Eq. [^. This suggests that the injection 
of twist via the emergence of current-carrying magnetic 
field is important for the creation of recurrent helical jets 
studied here. The amount of twist injection required be¬ 
tween succ essive jet-like episod es is consistent with the 
findings of Pariat et al. (2010), who carried out fully- 
compressible IVIHD simulations to model the formation 
of homologous helical jets. In their numerical experi¬ 
ments, twist injection is due purely to rotational mo¬ 
tion in the photosphere. However, MHD simulations of 


sphere (e.g 

. lArchontis & Hood 

|2013l IMoreno-Insertis & 

Galsgaard 

2013| 4'ang et al. 21 

)14| Lee et al. 2015) also 


of current-carrying field and Lorentz-force driven rota¬ 
tional motions) followed by emission of helical jets. It 
is likely both contributions are present for the observed 
helical jets studied here (which are found above a region 
of emerging flux). 

This work is an example of how complementary obser¬ 
vations from multiple observatories can be used in tan¬ 
dem with data-driven modeling to investigate the dy¬ 
namics of the solar atmosphere. Co-spatial and simul¬ 
taneous observations from IRIS, SDO and Hinode pro¬ 
vide evidence for the helical nature of the recurrent jets 
and reveal the magnetic environment of their source re¬ 
gion. The use of HMI vector magnetograms to perform 
data-driven simulations allowed us to investigate how 
processes such as flux emergence drive coronal evolu¬ 
tion. The data-driven simulations were carried out with a 
magnetofrictional model, which is able to capture how a 
magnetic configuration relaxes in response to the Lorentz 
force. However the model lacks substantial physics and 
is not suitable for answering questions related to how 
the stored magnetic energy is used to heat previously 
cool plasma to transition region and coronal tempera¬ 
tures, and how the plasma ejected along the jet is acceler¬ 
ated. Eor example, an important question regarding jets 
is whether the ascending material is directly accelerated 
by the Lorentz force in reconnected field lines, or whether 
it is due to chromospheric evaporation, sl ow mode waves, 


or up ward propagating shocks (e.g., see Takasao et al 
2013|). Some of these outstanding issues will likely be 


addressed in the near future by more detailed analyses of 
IRIS spectra of chromospheric and transition region lines. 
Going forward, an improvement over the present study 
would involv e data-driven, fully- compressible MHD sim¬ 
ulations (see Bonrdin et al.|2Q13 , for an example of data- 
driven MHD modeling applied to quasi-steady AR coro¬ 
nal loops). Measurements of the vertical gradients of the 
magnetic field from vector magnetograms at two heights 
(photosphere and chromosphere) will help retrieve the 
driving electric field, but advances in deriving consistent 
boundary conditions in terms of the appropriate mass. 
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momentum, and energy fluxes are also needed. The abil¬ 
ity to do so will allow us to strengthen the constraints 
imposed by observations (e.g. in terms of temperature 
and density diagnostics) on theory and to better inter¬ 
pret observations based on realistic physical models. 
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